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The notion of chiral symmetry for the conventional Dirac cone is generalized to include the tilted 
Dirac cones, where the generalized chiral operator turns out to be non-hermitian. It is shown that 
the generalized chiral symmetry generically protects the zero modes (n = Landau level) of the 
Dirac cone even when tilted. The present generalized symmetry is equivalent to the condition that 
the Dirac Hamiltonian is elliptic as a differential operator, which provides an explicit relevance to 
the index theorem. 

I. INTRODUCTION 

The chiral symmetry plays an important role in specifying some of the universality classes of the critical phenomena 
in disordered systems Systems are called chiral- symmetric when there exists an operator T that anti-commutes 
with the Hamiltonian, {T, H} = 0, with T^ = T and T 2 = 1. With this symmetry the energy eigenvalues appear 
always in pairs (E,—E), since if we have an eigenstate ?pE with an eigenvalue E, the state TipE is an eigenstate 
with an eigenvalue — E. The energy spectrum is therefore exactly particle- hole symmetric even when there exists a 
disorder as far as the disorder respects the chiral symmetry. In particular, the zero-energy state can be expressed as 
an eigenstate of T. 

For a massless Dirac fermions in two dimensions, usually with vertical and isotropic Dirac cones as in graphene^"— 
the effective Hamiltonian is expressed as Ho = vf(cf x t: x + cr y ir y ), where v x {y) is the Pauli matrix and 7r = p + eA 
denotes the dynamical momentum with the vector potential A and the electron charge — e. The fermi velocity is 
denoted by vp. For such a vertical Dirac cone, we have obviously {Ho,cr z } = and thus the system is chiral- 
symmetric with F = a z . The zero-energy Landau level (zero modes), which is essential to the anomalous quantum 
Hall effect for massless Dirac fermions in a magnetic field, then becomes an eigenstate of a z . The most remarkable 
property with these zero modes is their robustness against disorder in gauge degrees of freedom. The zero energy 
(n = 0) Landau level does not acquire any width due to such a disorder, while other Landau levels are broadened as 
usual, and this gives rise to an unconventional criticality for the quantum Hall transition at the n = Landau level. 7 
This robustness of zero modes for a vertical Dirac cone has been discussed in terms of the index theorem 8,9 or based 
on the explicit form of wave functions due to Aharonov and Casher j 10 ^ 1 It has been also demonstrated numerically 
that the chiral symmetry is also crucial to this robustness of zero modes of vertical Dirac cones.— 

Massless Dirac fermions in two dimensions appear not only in graphene but also in certain organic metals, where 
we encounter pairs of tilted Dirac cones iSr— The existence of zero modes and the Landau level structure have been 
established also for tilted Dirac cones. However, the effect of disorder, in particular the robustness of zero modes 
and the role of symmetry, has not been explored until recently^ In the present paper, with an explicit form of the 
eigenstates of the generalized chiral operator, we demonstrate how the chiral symmetry, which is broken for tilted 
cones, can be generalized to include tilted Dirac cones and clarify its relevance to the robustness of the zero modes of 
generic massless Dirac fermions. Relationships between the generalized chiral symmetry and the applicability of the 
index theorem is also elaborated. 

II. GENERALIZATION OF CHIRAL SYMMETRY 

To illustrate how the chiral symmetry is generalized to tilted Dirac cones, let us consider a general form of the 
Hamiltonian, 18 

H(rj) = -rjv F a 7T x + v F (a x 7r x + (T y 7r y ), 

for a two-dimensional massless Dirac fermion in a magnetic field, where the isotropic Dirac cone is tilted in the x 
direction for rj ^ with <jo being the unit matrix. The dynamical momentum 7r satisfies the commutation relation 
[ir x ,iTy\ = —iehB with B = rot A. The parameter rj determines the tilting of the Dirac cone. Note that for A = an 
equienergy contour of the Dirac cone is elliptic as long as \rj\ < 1 (while hyperbolic for |^| > 1; see FigJTJa)). The first 
term in the Hamiltonian destroys the chiral symmetry as a z H(r])a z = —H(—rj) ^ —H(rj). We can, however, define a 
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generalized chiral operator 7 as 

7 = \~ 1 {cr z - ir)a y ), = y/l -rj 2 , 

which arises naturally in the general framework of the eigenvalue problem for tilted Dirac cones ^ Although the 
generalized chiral operator 7 is not hermitian (7^ ^ 7), one can readily verify that 7 2 = 1 with eigenvalues ±1. The 
corresponding right-eigenvectors 7|±) = ±|±) are given explicitly as 

In the limit of the vertical cone (77 — >• 0), the generalized operator 7 reduces to the conventional chiral operator <j z . 
With the generalized chiral operator, we find that for |?7| < 1, 

7 f H(»7)7 = 

which we call the generalized chiral symmetry^ The generalized symmetry guarantees the identity (+\H(rj)\+) = 
(—\H(rj)\ — ) = 0, and plays an essential role for the robustness of zero modes as shown in the next section. Note that 
this symmetry holds irrespective of the details of the vector potential A. Disorder in gauge degrees of freedom (such 
as random magnetic fields) respects this symmetry. 

III. ROBUSTNESS OF ZERO MODES 

If the above generalized chiral symmetry is preserved, the Schrodinger equation H(rj)ijj = EiJj for the wave function 
ifj = + |— becomes 

/ \rj7Tx ~ iTTy \ / 0+ \ = E ( 1/Xrj T] / X^ \f ^\ 

\X v 7r x +i7Ty J J \ rj/Xr, J \ip J ' 

The zero (E = 0) modes are then given by the wave functions satisfying 18 

(X V 7T X - i7T y )^~ =0, ^ + = (1) 

or 

(XrqTTx + i7T y )^ = 0, 1/J- = 0. (2) 

The zero modes are thus the eigenstates of the generalized chiral operator 7 and have either "— " chirality with 
7*0 = — -0 (Eq.([TD) or chirality with ji/j = ip (Eq.©). It is to be recalled that for a vertical Dirac cone the zero 
modes are also the eigenstates of the chiral operator T = cr z . These equations for the zero modes hold even in the 
case where the gauge field is disordered. 

Following Aharonov and Casher^ we adopt the "Coulomb gauge" X^dxAx + X~ 1 d y A y = by assuming A = 
(-X^dyip, X v d x (p). Then Eq.© for ip~ and Eq.© for ip+ are reduced to 

[D ± T (e/fi)(D ± ^)]^ ± = 0, 

where D± = (dx ±idy) with R — (X, Y) — [xj y^\). The solutios are then given by ip^ = exp(±eip/h)f(Z±) 
with a polynomial f(Z±) in Z± = X ± iY . Since {d 2 x + d^)(p = B, we have <p(R) = j dR'G(R - R')B(R') where 
G(R) = (l/27r) log(ii) with R = VX 2 + Y 2 , which leads to the asymptotic form cp ->> (<I>/27r) logi? in the limit as 
R — » 00 with $ being the total magnetic flux in the system. We then have to chose the chirality for the wave function 
to be normalizable. For instance, when the total magnetic flux <I> is positive, the zero mode has to have "— " chirality 
to be normalizable. The number of square-normalizable wave function then becomes <I>/(/i/e), 10 that exactly coincides 
with the degeneracy of the Landau level. This implies the zero-modes exhaust the Landau level of a tilted Dirac cone 
so that its density of states is a generalized chiral symmetry protected delta-function in the presence of disorder. In this 
sense, the generalized chiral symmetry protects the zero mode of a generic massless Dirac fermions in two dimensions. 
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IV. GENERALIZED CHIRAL SYMMETRY AND INDEX THEOREM 

The above reasoning gives an explicit relationship between the index theorem and the generalized chiral symmetry. 
The index theorem, which gives the least number of zero modes, holds for an elliptic operator £ In the present case, 
the Hamiltonian H(rj) is elliptic as a differential operator if the matrix, 

is invertible for any (£ x , £ y ) G R 2 — (0, 0). The condition for the ellipticity thus becomes det H = — (1 — r] 2 )^ 2 — £ 2 ^ 0. 
This determinant is always negative (detS < 0) and never becomes zero for (£ X i€y) £ — (0,0), as l° n g as M < 1- 
The condition |^| < 1 is therefore exactly the condition for the ellipticity of the Hamiltonian H(rj) as well as that for 
the existence of the generalized chiral symmetry. In this sense, the generalized chiral symmetry is generically equivalent 
to the ellipticity of the Dirac-cone Hamiltonian with the same parameter space for their validity. Geometrically, the 
condition |?7| < 1 means that the Dirac cone is not tilted too much so that the cross section with a constant energy 
plane remains to be elliptic. 

V. NUMERICAL DEMONSTRATION 

We have also preformed numerical calculations based on a lattice model having a pair of tilted Dirac cones at E = 
as shown in Fig HJb) (Inset). The model is defined on the two-dimensional square lattice with a Hamiltonian having 
the nearest-neighbor (t) and the next-nearest neighbor (t f ) transfer integrals a o 18 i 2Q 

H = ^ ] ~t c r+e x c r H~ ( — 1) +2 ^<4, +e ^C r + t (c\.+ ex + ey C r + (? r J r&x - &y C r ) + H.C., 
r 

where r = (x, y) denotes the lattice point in units of the nearest-neighbor distance a, and e x (e y ) the unit vector in the 
x(y) direction. The magnetic field is taken into account by the Peierls phases 0{r) as t(t') — »• t{t') exp(— 2iti9(r)) 1 so 
that the their sum along the closed loop is equal to the enclosed magnetic flux in units of the flux quantum h/e. To see 
how the Landau levels are broadened by disorder in the present model, we assume that the magnetic flux <p(r) enclosed 
by the square located at r has a random component 5<j)(r) in addition to the uniform component </>. The random 
components have a gaussian distribution and are correlated in space as {5<j){ri)8(j){r2)) = {8(j) 2 ) exp(— \r\ — T2 1 2 /4?7^). 
The disorder in magnetic fluxes should appear, for large 77^, as the disorder in gauge degrees of freedom, and therefore 
should respect the generalized chiral symmetry of the effective Hamiltonian at the Dirac points. As shown in FigOJb), 
we actually see an anomalously sharp n = Landau level when the correlation length of disorder becomes larger 
than the nearest-neighbor distance a, while other Landau levels (n — ±1,±2,±3) are broadened as usual. Together 
with the result 18 for a different inclination of the cone, this anomaly at the n = Landau level is likely to exist 
universally, irrespective of the magnitude of the tilting angle of the cone. This anomalously sharp n = Landau level 
suggests that the energy levels are exactly degenerated at E = 0. The present numerical results are thus consistent 
with the robustness of zero modes of a single tilted Dirac cone protected by the generalized chiral symmetry. 

VI. CONCLUSIONS 

We have shown explicitly that the notion of the chiral symmetry can be generalized to generic tilted Dirac cones, 
where the generalized chiral operator has to be non-hermitian. It has been also demonstrated analytically that the 
generalized chiral symmetry indeed protects the zero modes of the system by extending the argument by Aharonov 
and Casher*i£ The resulting anomalously sharp n = Landau level has been also confirmed numerically based on the 
lattice model. The existence of the generalized chiral symmetry coincides with the ellipticity of the Hamiltonian as 
a differential operator, which is nothing but the geometrical condition that the Dirac cone is not tilted too much so 
that the cross section with a constant energy plane is an ellipse. 
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FIG. 1: (a) Schematic figures of the Dirac cones with 77 = 0, \rj\ < 1 and \rj\ > 1. (b) Density of states for the lattice model 
in a magnetic field with various values of the correlation length 77^ of the random component of the magnetic flux. Here, the 
parameters are assumed to be t' /t = 0.25, (j) = 0.01(/i/e), and yj ' {8(j) 2 ) = 0.0029(/i/e). The system -size is 20a by 20a and the 
average over 10 4 samples is made. Inset: Energy dispersion E(k)/t in the absence of a magnetic field, where k\ — k ■ (e x + e y ) 
and &2 = k • (e x — e y ). 
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